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evaluated using the solar radiation data collected by the meteorological NRG weather station installed 
in the Center of Research and Technologies of Energy (CRTEn) in the Borj-Cedria area. The collected 
measurements during the last three years (2008, 2009 and 2010) were based on 10 min time step. 
Keywords: These data have allowed us to evaluate the global solar flux, the sun duration, the yearly and the 
Gulf of Tunis seasonal frequency distribution of the global solar radiation. Moreover, a conventional model has been 
Global solar Radiation used to estimate the hourly solar radiation on a horizontal plane and it has been validated by 
Model comparison $ f 5 eels 3 
Solar energy potential experimental measurements in specific days. The results show that the global solar radiation predicted 
by the conventional model has a good agreement with the experimental data during the clear sky 
conditions with a mean absolute percentage error (MAPE) of 4.1%. However, the limitation of the 
conventional model appears under the cloudy sky weather which is proved by the highest value of 
relative error percentage reaching 14.26% occurred during the autumnal equinox day. 
© 2012 Elsevier Ltd. All rights reserved. 
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1. Introduction 


No doubt that all kind of energies took origin from the sun. 
This huge furnace, with a diameter of 1.39 x 10° m, emits a total 
energy output of 3.8 x 102° MW, however only a tiny fraction of 
1.7 x 10'4kW is intercepted by the earth. So that since the 
prehistoric times, man has realized that a good use of solar 
energy is in his benefit. Therefore, solar energy has been har- 
nessed by mankind [1]. Indeed, the Greek physician Archimedes 
(212 BC) burned the Roman fleet by means of concave metallic 
mirror in the form of hundreds of polished shields; all reflecting 
on the same ship [1,2]. Then, Kircher (1601-1680) worked out 
some experiments to validate the story of Archimedes. Yet, no 
report of his findings survived [3]. 

Despite of the serious difficulties to apply the solar energy 
concentration, the use of concentrating collectors was the first 
application developed in this field. The construction of the first solar 
furnaces, capable of melting iron, copper and other metals, was 
recorded from the French scientist Lavoisier during the 18th 
century [1]. 

The furnaces were in use throughout Europe and the Middle 
East and attained the remarkable temperature of 1750 °C. The 
19th century was marked by two important events; the first one 
was the discovery of the photovoltaic effect by Becquerel as a new 
form of solar energy exploitation. The second was related to the 
several attempts to convert solar energy into other forms based 
upon the generation of low pressure steam to operate steam 
engines. In fact, Monchot has built several solar-powered steam 
engines between the years 1864 and 1878 [3]. 

As the lack of water was always a problem to humanity, 
researchers attempted to harness solar energy by developing new 
equipment suitable for the desalination of sea-water. Conse- 
quently, solar distillation has been in practice for a long time. 
At the same period, the first water distillation plant was built at 
Las Salinas, Chile in 1874 [3,4], using solar stills that evaporate 
salty water. Besides, in Algeria Mouchot designed and set a 
truncated cone reflector in1875. Another furnace was set up in 
Algeria. 

After that, Eneas installed a focusing collector to operate a water 
pumping apparatus at California farm in 1901 by concentrating the 
sun’s rays at a focal point where the boiler was located. In 1904, a 
Portuguese priest, constructed a large solar furnace which appeared 
quite modern in structure, being a large, off-axis, parabolic horn 
collector [3]. 

In 1913, the world’s largest pumping plant was implemented in 
Meadi, Egypt by Shuman, in collaboration with Boys using long 
parabolic cylinders to focal point sunlight into a long absorbing tube. 
In 1928, Pasteur has reported a novel way to use solar concentrators 
in solar distillation by focusing solar rays into a copper boiler 
containing water. The steam generated from the boiler was piped 
to a conventional water cooled condenser in which distilled water 
was accumulated [4]. 

Otherwise, the house heating and hot water appeared in the mid 
1930s and recognized a real expansion in the last half of the 40s all 
over the world. The main solar water heaters were manufactured in 


many countries, for instance, the thermosyphon types which consist 
to a two flat-plate solar collectors having an absorber and storage 
tank. 

On the other hand, the first practical application of solar cells 
was in space in the 1960s when scientists discovered other 
photovoltaic materials such as gallium arsenide (GaAS) which 
increased the cost of PV installation. 

The global installed capacity of photovoltaics at the end of 
2002 was near 2 GWp [5] Photovoltaic (PV) cells are made of 
various semiconductors, basically the silicon (Si), compounds of 
cadmium sulfide (Cds), cuprous sulfide (Cu2S), and GaAs which 
are materials that are only moderately good conductors of 
electricity. Thirty years ago, PV modules began to be for large- 
scale commercial applications with an efficiency below 10%. 
Nowadays, their efficiency is about 20% and estimated to reach 
25% in 2020. 


2. Recent studies 


Knowledge of the availability of the global solar flux at the 
earth’s surface has a great importance in solving many scientific 
problems and harnessing the practical utilization of solar energy. 
This available quantity of energy changes within geographic place 
because it depends on climatic and atmospheric conditions. 

It is known that global solar radiation is poorly sampled in 
weather station networks. For this reason, we find in the literature, 
recently, many authors who attempted to develop several models to 
estimate the global solar radiation and asses their corresponding 
solar potential [6-23]. We can find also software packages allowing 
the solar flux density estimation of several sites. 

Coskum et al. [24] have modified the concept of the probability 
density frequency usually used to analyze the wind and the 
outdoor temperature distributions in order to estimate the solar 
radiation distribution and therefore to investigate and design 
better the solar energy systems. They have also conducted their 
study with global solar irradiation data of many years recorded by 
the Turkish State Meteorological Service. 

Ryder et al. [25] have used the solar irradiance measurements 
from a new high density urban network in London. They have 
measured annual averages and demonstrated that central London 
receives 30 + 10 W m~? less solar irradiance than outer London at 
midday, equivalent to 9 + 3% less than the London average. They 
have obtained these measurements basing on a new technique 
referred to the ‘Langley flux gradients’ that infer aerosol column 
concentrations over a clear periods of 3 h and this technique has 
been developed and applied to three case studies. 

The survey of Ramachandra et al. [26] has focused on the 
assessment of the Indian potential resource with variability 
derived from high resolution satellite data. They have also 
presented a techno-economic analysis of the solar power technol- 
ogies and a prospective for a minimal utilization of the available 
land. 

Sorapipatana [27] has adopted a satellite technique to assess 
solar energy potential in Kampuchea. In his work, he has estimated 
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the solar irradiation potential at an interval of half a degree grid. The 
seasonal variations of mean daily solar irradiation in Kampuchea 
were measured during two Asian winter and summer monsoon 
seasons. 

Molero et al. [28] have developed a transient 3-D mathema- 
tical model for solar flat plate collectors. Their model is based on 
setting mass and energy balances on finite volumes and allows 
the comparison of different configurations. This model is a useful 
tool to improve the design of plate solar collectors and to compare 
different configurations. 

Hocaoglu [29] proposed a solar radiation modeling approach 
using the hidden Markov models. This approach mean of hourly 
measured ambient temperature values are considered as observa- 
tions of the model, whereas mean of hourly solar radiation values 
are considered as the hidden events, which constitute the out- 
comes of the proposed mathematical model. The proposed model 
exploited the inherent cross dependency between daily ambient 
temperature and solar radiation in a form of hidden Markov 
model (HMM) using the optimization technique called the Viterbi 
decoding algorithm. 

Ozgoren et al. [30] developed an artificial neural network 
(ANN) model based on multi-nonlinear regression (MNLR) 
method for forecasting the monthly mean daily sum global solar 
radiation at any place of Turkey. This aim was based on the 
meteorological data of 31 stations spread over Turkey along the 
years 2000-2006. The developed model was very useful to 
calculate quickly and satisfactorily the monthly global solar 
radiation at any place of Turkey. 

Maghrabi [31] established a simple model to estimate the 
monthly mean global solar radiation on a horizontal surface in 
Tabouk (Saudi Arabia) using 9 years of solar radiation data. The 
model correlated the global solar radiation with five meteorolo- 
gical parameters. These parameters were the perceptible water 
vapor, air temperature, relative humidity, atmospheric pressure, 
and the mean monthly daily fraction of possible sunshine hours. 
The obtained model showed its high accuracy for the monthly 
global solar radiation estimation in Tabouk and other places 
having similar conditions. 

Mao and Li [32] evaluated the monthly mean daily solar 
radiation in Chongqing basing on the observed meteorological 
variables including daily maximum and minimum temperatures, 
daily mean dew point temperature, fog and rainfall analyzed from 
1986 to 2000. Moreover, a newly developed model that included 
all selected variables proved the best method with a root mean 
squared error (RMSE). In fact, 12 models were developed however 
they performed differently in different seasons. 

Moradi [33] developed a new and automatic method for 
controlling the quality of daily global solar radiation using 
sunshine duration hours. In fact, the proposed method was based 
on a three levels of tests: which were firstly the daily global solar 
radiation has to be compared against the daily extraterrestrial 
radiation that is received on a horizontal surface second, secondly 
the daily global solar radiation should only exceed by a small 
amount of the daily clear sky irradiation that was observed under 
highly transparent clear skies and thirdly the method used a 
series of persistence checks that utilized the relation between 
daily global solar radiation and relative sunshine duration hours. 
The method proved its ability e to identify the systematic and 
non-systematic errors under three different climates including 
semi-arid, coastal humid and very arid climates. 

Voyant et al. [34] presented an application of Artificial Neural 
Networks (ANN) to predict the daily solar radiation. The authors 
proposed a contribution study of exogenous meteorological data 
compared with different forecasting methods. This study showed 
that the use of exogenous data generates a RMSE between 0.5% 
and 1% for two stations during 2006 and 2007 (Corsica Island, 


France). Besides, the prediction results were also relevant for the 
concrete case of a tilted PV wall. In fact, the addition of 
endogenous and exogenous data allowed a 1% decrease of the 
RMSE over 6 months-cloudy period for the power production. 
While the use of exogenous data showed an interest in winter, 
endogenous data as inputs on a preprocessed ANN appeared 
sufficient in summer. 

Korachagaon and Bapat [35] developed a general formula to 
estimate the solar radiation with the local site parameters such as 
air temperature, relative humidity, wind, moisture and few 
derived parameters independent variables. The developed models 
were validated with remaining 665 data sites. Authors concluded 
that two candidate models should be proposed because they were 
are capable to cover 50% of the land area on earth surface 
between latitude + 30 with an estimation error (RMSE) limiting 
to 15%. Thus, it was predicted that the proposed models can be 
used to estimate the monthly average daily global solar radiation 
in area where the radiation data is missing or not available. 

Si-qing et al. [36] made a study on the medium-term forecast 
of solar 10.7 cm radio flux (F10.7). By comparing their forecast 
results in the period from 21th September 2005 to 7th June 2007, 
they have been demonstrated that the accuracy of the autore- 
gressive forecast method is equivalent to that of the forecast 
made by the American Air Force. 

In this study, an overview of research works of all over the 
global solar radiation models is presented. The mentioned models 
can be considered as the representative of worldwide average. Its 
details can be found therein. Moreover, we present a classic 
theoretical model of the solar flux density of sun in the site of 
the gulf of Tunis in Tunisia. This model is used to estimate the 
direct and the diffuse radiation. 

We proceed to the validation of this model by experimental 
measurements which are conducted by means of a meteorological 
station during the period of 2008-2010 with a 10-min time step 
by the NRG weather station. The treatment of about 158,000 
observations permit us to evaluate the global solar flux in the site 
of the gulf of Tunis. 


3. Site description and experimental design 


Tunisia is located in North Africa. It is bordered by Algeria to 
the west and Libya to the southeast. It is the northern most 
country on the African continent, and the smallest of the nations 
situated along the atlas mountain range (Fig. 1). The country is 
divided in two regions, the well-watered north and the semi-arid 
south. An abrupt southern turn of its shoreline gives Tunisia two 
faces on the Mediterranean Sea with 1298 km coastline. Fig. 1 
shows also that Tunisia is situated in the most insolated regions 
of the globe for which the annual average global solar radiation 
exceed 2000 kWh/m7?/year. Fig. 2 shows the geographic localiza- 
tion of the Borj-cedria site in the gulf of Tunis. 

The experimental design is a NRG station installed in 
36°43'04’N latitude and 10°25’41”E longitude (Fig. 3). It provides 
the measurement of wind speed, wind direction, ambient tem- 
perature and the global solar flux in the gulf of Tunis [37]. This 
station is equipped with an acquisition system that record, every 
10 min the average, the max, the min and the standard deviation 
values for each sensor. 


4. Theoretical models of solar flux estimation: an overview 


In the literature, several empirical models have been developed 
to estimate global solar radiation which depends on different 
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Fig. 1. Geographic position of Tunisia. 


factors such as extraterrestrial radiation, sunshine hours, mean 
and maximum temperature, relative humidity, number of rainy 
days, altitude, latitude, the turbidity of the atmosphere and so 
forth [38-61]. In general, most of the widespread literature 
available on evaluating the surface solar flux can be categorized 
into two broad groups: monthly-averaged daily irradiation and 
monthly-averaged hourly irradiation. 


4.1. Monthly-averaged daily radiation 


Measured monthly-averaged values of daily irradiation H are 
an excellent basis of information and offer the starting point of 
many computation schemes. Angstrom (1924) proposed a corre- 
lation relating H and the monthly average of the instrument - 
recorded daily time fraction of bright sunshine S as follows [62]: 


H; =%+0-0S (1) 
where Ho and «& are the monthly average daily extraterrestrial 
radiation and location-based empirical constant (#=0.25 for Stock- 
holm) respectively. Further efforts in the literature suggested 
additional empirical modifications to the Angstrom equation are 
illustrated in the following section. 


4.1.1. Angstrom—Prescott-Page model 

The Angstrom-Prescott-Page model is the most commonly 
used model as given by [63]: 
H S 
Ho =a+ a (2) 
where H, Ho, S, So, a and b are the monthly average daily global 
radiation, the monthly average daily extraterrestrial radiation, the 
day length, the maximum possible sunshine duration, and two 
empirical coefficients respectively. 


Fig. 2. Geographic localization of Borj Cedria site. 


The monthly average daily extraterrestrial solar radiation 
on a horizontal surface can be given using the following 
equation [64,65]: 


24 ‘ T ; er 
Ho = g on [cos (cos Osin(Wws) + 730 wssin ọsin J (3) 
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Fig. 3. NRG weather station. 


where Gon is normal extraterrestrial radiation, ~ is the latitude of 
the site, ô is the solar declination and w, is the mean sunrise hour 
angle for the given month. The normal extraterrestrial radiation, 
the solar declination and the sunrise hour angle can be computed 
by the following equations: 

360 2) 


Gon = Gsc (1 +0.033cos ——_— 


365 4 


where G,c is solar constant which is defined as the energy from 
the sun per unit time received on a unit area of surface 
perpendicular to the direction of propagation of the radiation at 
mean earth-sun distance outside the atmosphere and n is the 
ordinal date. n=1 at January 1st and n=365 at December 31st. 


ae 284+n 
ô =23.45°sin (360 oo | (5) 
ws = arcos(—tan otan ô) (6) 


For a given month, the maximum possible sunshine duration 
So can be estimated as follows [66]: 


2 
So = 5 Ws (7) 


4.1.2. Tiris et al. model 

Tiris et al. proposed the a and b empirical coefficients of the 
Angstrom-Prescott-Page model for Turkey, in general, as follows: 
H S 


Ho lee (8) 


4.1.3. Page model 

Page has given the coefficients of the Angstrom-Prescott-Page 
model, which is believed to be adopted all over the world, as the 
following [67]: 
H S 


H; 7023+048 9) 


4.1.4. Bahel et al. model 
Bahel et al. proposed the following relationship [68]: 
S 


H 
Ho Soe (10) 
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4.1.5. Louche et al. model 

Louche et al. suggested the following model to expect the 
global solar radiation [69]: 
H S 


H = 0206+0.546 > (11) 


4.1.6. Benson et al. model 
Benson et al. presented two different formulations for two 
periods of a year as the following [70]: 


H S 

m 7018+060 a2) 
for the periods during January-March and October-December. 
H S 

Ho SORAR Uae (13) 


for the period between April and September. 


4.1.7. Monthly specific Dogniaux and Lemoine model 
Dogniaux and Lemoine proposed the following specific 
monthly correlations [71]: 


January “+ = (0.003016 +0.34507) + (0.004958 +.0.34572) Å 
0 (0) 
a4 


February = (—0.003016 + 0.34507) + (0.004958 + 0.34572) $ 
0 0 


(15) 
H : S 
March 7- = (—0.003038 + 0.36690) + (0.004666 +0.36377) ©- 
0 ‘0 
(16) 
-H S 
April 7 = (-0.003345 + 0.38557) + (0.004566 + 0.35802) = 
0 0 
(17) 


May jr = (0.002455 +0.35057) + (0.004855 +0.33550) = (18) 


June $ = (0.003275 +-0.39890) + (0.005785 +-0.27292) > (19) 
0 10) 


July f = (0.003695 +-0.41234) + (0.005685 +0.27004) > (20) 
0 0 
H S 
August z- = (—0.002695+ 0.36244) + (0.004125 +0.33162) ©- 
0 0 
(21) 


September F- = (0.003386 +0.39467)+ (0.005645 +0.27125) > 
0 0 


October -+ = (0.003716 +0.36213)+ (0.005045 +0.31790) Å 
0 10) 
(23) 


November- = (—0.003506 + 0.36680) + (0.005238 + 0.31467) > 
0 0 
(24) 


December jn = (—0.003506 + 0.36262) + (0.005596 + 0.30675) 7 
0 0 


(25) 


4.1.8. Specific Monthly-Rietveld model 
Soler used the Rietveld’s model to 100 European stations and 
suggested the following specific monthly relationships [72]: 


H S 
January z =0.18 +0.66 5 (26) 
February -L = 0.20+0.60— (27) 
Ho So 
H S 
March yy = 0.224038 5 (28) 
Ht S 
April z = 0.20 + 0.62 So (29) 
H S 
May i = 0.24 + 0.52 5 (30) 
H S 
June z =O AHU (31) 
H S 
July =0.23+0.53 5 (32) 
H S 
Augusty = 02240555 (33) 
H S 
September — = 0.20+ 0.59 — (34) 
Ho So 
October cS = 0.19+0.60— (35) 
Ho So 
November- =0.1 740.662 (36) 
Ho So 
December =0.18+0.65 > (37) 
Ho So 


4.1.9. Gopinathan model 
This model assumed to use the following Angstrom-Prescott- 
Page coefficients [73]: 


a=0.265+0.07Z—0.135 (=) (38) 
(0) 
S 
b = 0.401 +0.108Z—0.325 (=) (39) 
0 


where Z is the altitude of the considered site in kilometers [64]. 


4.1.10. Zabara model 

Zabara suggested monthly values to the coefficients a and b of 
the Angstrom-Prescott-Page model with monthly relative sun- 
shine duration (s/So) as a third order function and expressed a and 
b as follows [74]: 


S 542 s\3 
a—0.395—1.247 + 2.680 1.647 (40) 

So So So 

S Sy 53 
b=0.395+1.384 3.249 +2.055 (41) 

So So So 


4.1.11. Kilic and Ozturk model 

Kilic and Ozturk provided empirical values of the a and b 
parameters of the Angstrom-—Prescott-Page model for Turkey as 
follows [75]: 


a= 0.103+0.000017Z + 0.198cos(p—6) (42) 
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b=0.533+0.165cos(p—6) (43) 


4.1.12. Gariepy’s model 

Gariepy has expressed empirically the coefficients a and b of 
the Angstrom-Prescott-Page model in terms with mean air 
temperature T and the amount of precipitation (P) [76]: 


a= 0.3791—0.0041T—0.0176P (44) 


b=0.4810+0.0043T + 0.0097P (45) 


4.1.13. Almorox and Hontoria model 
Almorox and Hontoria suggested an exponential model as 
follows [77]: 


H S 
H = a+bexp (=) (46) 


4.1.14. Ampratwum and Dorvlo model 
Ampratwum and Dorvlo proposed the following model [78]: 


H S 
i= a+ blog (=) (47) 


4.1.15. Ogelman et al. model 
Ogelman et al. have correlated (H/Ho) with (S/So) in the form of 
a second order polynomial equation [79]: 


H S Si 

(a) Hy =a+b(=) (=) (48) 
H S S\? 

(b) i= 0.195+0.676 (=) —0.142 (=) (49) 


where a, b and c are experimental coefficients. 
4.1.16. Samuel model 


The ratio of global to extraterrestrial radiation was expressed 
by a function of the ratio of sunshine duration as follows [80]: 


H S Ss\? sy 
@) Ho aA (=) (5) ee G Ho 


H S S\? s\3 
() p =0.1442.52( =.) 3.71(¢-) + 2.24 (=) (51) 


where a, b and c are empirical coefficients. 


4.1.17. Allen model 

Allen developed a self adjusting model which is a function of 
the daily extraterrestrial radiation, mean monthly maximum and 
minimum temperature [81]: 
H 
Ho 


= A(T max Ta +b (52) 


4.2. Monthly-averaged hourly radiation 


The second group of methods deals with the estimation of the 
monthly-averaged hourly solar-radiation. One of the most used 
models among this category is the ASHRAE clear-sky model [82]. 
In this model, the direct normal irradiation is estimated by means 
of a simple equation including two constants A and B while the 
diffuse radiation is given as a fraction C of the direct normal 
component. 

The hourly global solar radiation on horizontal plane Q, is 
composed by a beam (direct) radiation Q,, and a diffuse sky 
radiation Qan. 


The beam radiation Qp, is given by the following formula [82]: 
Qbh = Qn x cos0z (53) 


where Q, and 0, are the normal direct radiation and 0; is the 
zenith angle respectively. 
The diffuse sky radiation Qan can be expressed as yields: 


Qan = Qn x C (54) 


where C is the diffuse-sky factor. 
Therefore, the global radiation on horizontal surfaces is: 


Qan = Qn x cos0z+ Qn x C (55) 


We note that direct normal radiation Q, is calculated as the 
following: 


Qn =A x exp(—B/sinas) (56) 


where A is the apparent solar-radiation constant, B is the atmo- 
spheric extinction coefficient, and &s is the solar altitude angle 
above the horizontal (the complement of the zenith angle). 

ASHRAE [82] gave 12 sets of these three constants for each 
month of the year. The different sets of these constants are related 
to the monthly variation of the amount of precipitate water 
during the year. Hence, the twelve sets of ASHRAE’s coefficients 
reflect the annual variation of the absolute atmospheric humidity 
[83]. In fact, the constant A is related to the solar constant while 
the variations of the constants B and C indicate a variation in 
turbidity as well. In order to account for regional variations of 
humidity and turbidity, ASHRAE published maps for a parameter 
called “clearness index”, for both summer and winter, for differ- 
ent regions in the USA. This factor changes the values of amount 
radiation estimated by the model. Unfortunately, the unavail- 
ability of these factors for other regions of the world prohibited 
the use of this model for these regions. For this reason, several 
research works are focused to develop adjustment factors to the 
ASHRAE clear-sky model for many countries over the entire 
world. For instance, Powell [84] examined The ASHRAE model 
basing on data collected at 31 NOAA (National Oceanographic and 
Atmospheric Administration, USA) monitoring stations. The 
results established the universal validity of the model in predicting 
the solar radiation under cloudless conditions. Powell modified the 
basic ASHRAE model using altitude corrected optical air mass 
instead of seasonal clearness indexes. The author concluded that 
his modifications made the model generally more truthful. Machler 
and Iqbal [85] used the above modifications made by Powell to the 
ASHRAE model and adjusted the constants A, B, and C. Also, they 
introduced a correction humidity factor to account the variability of 
the water-vapor absorption. 

Galanis and Chatigny [86] criticized some formulations of the 
ASHRAE model. The authors proposed to include the clearness 
index in the expressions of the direct and more importantly in the 
diffuse radiation under cloudless conditions. Moreover, they 
suggested to re-write expressions for cloudy-sky conditions in 
order reduce to the cloudless formulation for zero cloud cover. 

Maxwell [87] developed the METSTAT solar radiation model 
(called METSTAT model) basing on data collected by 29 US 
National Weather Service sites. The model estimated the hourly 
values of direct normal, diffuse, and global solar radiation. The 
model requires some inputs such as the aerosol optical depth, the 
surface albedo, the snow depth, the atmospheric pressure, the 
precipitable water vapor and so forth. Even though, this model 
seems to be comprehensive, its application is restricted to regions 
where the above input data are available. 

Rigollier et al. [88] developed another clear-sky model in the 
framework of the new digital European Solar Radiation Atlas 
(ESRA). The factors in the model have been empirically adjusted 
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by fitting techniques basing on hourly measurements over several 
years for a number of European regions. 

Gueymard [89] suggested two models to calculate the 
monthly-average hourly global irradiation distribution from its 
monthly-averaged daily counterpart. He claimed that a quadratic 
expression in the sine of the solar elevation angle fitted accurately 
the data at all considered locations. The author included other 
parameters in his proposed models such as the mean monthly 
clearness index, the average day-length, and daily average solar 
elevation. He used a large data set from 135 stations covering 
diverse geographic locations and he confirmed the performance of 
those models presenting better relative precision compared to 
other previous models. 

Yang and Koike [90] suggested a numerical model to predict 
the hourly-mean global solar-radiation taking account the upper- 
air humidity. The authors defined a sky-clearness indicator as a 
function of the relative humidity profiles in the upper atmo- 
sphere; then they used this indicator to relate global solar 
radiation under cloudy skies to that under clear skies. 


4.3. Conventional way of estimating hourly solar radiation 


The following calculations of the global solar flux on a 
horizontal plane are based on a conventional model [91] which 
assumes that the amount of direct radiation, J, intercepted by a 
perpendicular surface to the sun rays which are expressed as 
follows: 


I= Gsc x exp| — 


TL (57) 


0.94+.9.4sina, 
where G,,, TL and as are the solar constant, the Linke’s trouble 


factor and the solar altitude respectively. The Linke’s trouble 
factor TL is expressed as the following: 


TL=2.4+14.6f +0.4(1 +2f)In(pp) (58) 


where f and pp are the atmospheric trouble and the partial 
pressure respectively. 

The value of f depends on atmospheric conditions of the site 
in question; 6=0.05 in rural zone, /=0.1 in urban zone, /=0.2 in 
industrial park or polluted. 

The partial pressure of steam p, (expressed in mm Hg) can be 
calculated by the following equation: 


760 
Pp = 701325 


where Pat and Par are the atmospheric pressure (in Pascal) and the 
pressure of the dry air (equal to 101222.105 Pa) respectively. 

The solar altitude «s, which is the angle between the horizon- 
tal and the line to the sun, can be given using the following 
relation: 


(Par—Par) (59) 


sing, = sin osin 6+Ccos cos 6cos w (60) 


where ¢ is the latitude of the site which is the angular location 
north or south of the equator, north positive — 90° < p < +. 90°. 
The hour angle œ is the solar angular displacement east or 
west of the local meridian due to the rotation of the earth (15° per 
hour). The hour angle is negative in the morning and positive 
afternoon. Therefore, the hour angle can be written as follows: 


@ = (ts—12)15° (61) 


where ts is the solar time. It is known that converting civil time, 
tc, to solar time ts, requires that two corrections should be 
applied; the first one is the correction for the difference in 
longitude 2 between the observer’s meridian and the meridian 
on which the civil time is based, the second correction is from 
the equation of time ET. Solar time and equation of time are 


expressed by Eqs. (62) and (63): 
À 
ts= tc+ J57 Zc+ ET (62) 


where Zc is the time zone east of GMT. 


ET = — (0.0002—0.47497cos(B) + 3.2265cos(2B) 
+0.0903cos(3B) + 7.3509sin(B) + 9.392sin(2B) 
+ 3.3361sin(3B)) (63) 


where B is the mean anomaly of the earth, called also the 
fractional year, which is a parameter relating position and time 
for the earth moving in a Kepler orbit. It is what the true anomaly 
would be if the earth moved with constant speed along a perfectly 
circular orbit around the sun in the same time. 

According to this model, the beam solar radiation component 
Gp received by a horizontal surface is given by the following 
equation: 


Gp = Gsc x exp Li | x sing, = Í x sings (64) 


~ 0.9+9.4sine, 
The diffuse solar radiation component Gy is calculated by the 
following formula: 


Ga = 54.8 \/sina.[TL—0.5—sines] (65) 


The global solar radiation on horizontal surface, Gr, is the sum 
of the diffuse radiation, Gg and the beam radiation. Therefore, the 
incident global radiation on horizontal surface is given by the 
following expression: 


TL 
Gr = Ga + Gp = Gsc x EXP -SSF asinw 
+54.8 \/sina,[TL—0.5—sings] (66) 


5. Experimental results 


In this section of study we are interested to the treatment of 
data registered by the NRG weather station installed on the site. 
This meteorological station is equipped by very precise and 
reliable instruments allowing the measure of the global stream 
of solar energy on a horizontal plane by way of a pyranometer 
(class A), the statement of ambient temperature and the sun 
covering rate. The programmed time step of the measurements is 


Table 1 
Ambient temperature characteristics. 


Ambient Temperature (°C) 


2008 2009 2010 
January 11.88 11.81 12.71 
February 11.96 11 12.81 
March 13.81 13.4 8.33 
Apri 16.95 15.34 17.16 
May 19.89 19.72 19.23 
June 23.29 23.29 22.28 
July 26.62 26.39 26.25 
August 26.66 27.84 26.28 
September 24.32 24.14 23.76 
October 21.38 19.96 20.72 
November 16.20 17:37 16.60 
December 11.87 17.19 13 

Minimum Ambient Temperature (°C) 

1.517 1.529 1.519 

Maximum Ambient Temperature (°C) 

40.37 41.33 41.17 

Average Ambient Temperature (°C) 

18.70 18.94 21.97 
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10 min. We note that solar radiation and sunshine duration has 
been abundantly measured in the country. A few attempts to 
develop knowledge about the variation of global solar radiation in 
Tunisia have been found in the literature. Thus, the aim of this 
paper is to assess the potential of solar energy resource in the gulf 
of Tunis. The statically measurements allow us to evaluate the 
average sun duration and the average solar flux during the years 
2008, 2009 and 2010. Also, the accuracy degree of the conven- 
tional model to the experimental data is investigated. 


5.1. Ambient temperature and global solar flux 


Table 1 recapitulates the characteristics of the ambient tem- 
perature during the last three years. These measurements 
(157,860 samples) are provided by the NRG weather station with 
a step of 10 min. The variation of the ambient temperature during 
the three years is shown in Fig. 4. It is observed in the first six 
months of each year that the ambient temperature has the trend 
to increase from the lowest temperatures in December surround- 
ing 10°C to the highest temperatures in the summer months 
surrounding 30 °C. However, in the second half of the year the 
ambient temperature decreases till the end of the year. It is worth 
noting that the ambient temperature has nearly the same trend of 
variation during the three years. One can observe also that during 
the three years the average ambient temperature increases from 
18.70 °C in 2008 to 21.97 °C in 2010. 

Table 2 illustrates the variation of the global solar radiation on 
a horizontal surface during the last three years. It is observed that 
the peak of average monthly radiation in the gulf of Tunis is 
registered in July with the value of 230 kW/m? while the least is 
December with the value of 63 kW/m?. The total annual global 
solar-radiation received in the Gulf of Tunis on a horizontal 
surface is 1681, 1763 and 1704 kWh/m? during 2008, 2009 and 
2010 respectively. It is noticeable that the sunshine durations are 
important with the values of 3337, 3505 and 3415 h during 2008, 
2009 and 2010 respectively. One can observe also that that more 
than 38% of this total global solar flux is contributed in the 
summer months. 


Ambient Temperature (°C) 


0 5000 10000 


Fig. 5 shows the monthly variation of global solar flux as 
hourly values, one can deduce obviously that the global solar flux 
has the same trend of variation as the ambient temperature 
presented in Fig. 4. Also, we find the highest global solar radiation 
values occurs in the summer months (June, July and August), 
while the lowest values lie in the winter months (December, 
January and February) and the global solar radiation values are 
seen to be equal through the autumn and spring months. This is 
because in the northern hemisphere, the sun’s elevation angle 
during winter (November, December, and January) is lower than 
during the other seasons [92]. It is explained by the fact that the 
mean distance between the Earth and the Sun in winter is less 
than that in spring (or autumn) by about 3% due to the eccen- 
tricity of the Earth’s orbit, resulting in a stronger flux of radiation. 
One can deduce that the variation of the global solar flux has the 
inverse trend of the zenith angle variation of the sun during the 
year. In the gulf of Tunis and in other mid-latitude countries north 
of the equator, the sun’s daily trip is an arc across the southern 
sky. The sun’s greatest height above the horizon occurs at noon, 
and how high the sun then gets depends on the season of the 
year, it is the highest in mid-summer (June), which correspond 
to the lowest zenith angle value, the lowest in mid-winter 
(December) that correspond to the highest zenith angle value. 


Table 2 
Global solar radiation characteristics. 


Period 2008 2009 2010 Average/month 
January 81 62 69 71 
February 98 84 95 92 
March 133 145 127 135 
April 179 155 157 163 
May 187 230 187 201 
June 220 244 217 227 
July 223 238 228 230 
August 189 203 202 198 
September 122 142 138 134 
October 107 109 110 109 
November 81 86 76 81 
December 61 65 64 63 
Average 140,109 146,992 139,233 142,111 


15000 20000 25000 


Obsercavtions 


Fig. 4. Ambient temperature variation. 
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Fig. 5. Monthly variation of the global solar radiation. 
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Fig. 6. Yearly frequency distribution of global solar radiation. 


Hence, this trend of variation can be related to how high in the 
sky the sun appears. 


5.2. Yearly and seasonal probability distribution of global solar 
radiation 


Fig. 6 shows the yearly probability distribution or the fraction 
of time, based on the averaged hourly measurement of the global 
solar radiation of the three last years where they are within the 
interval given by the width of the columns (here 50 W/m? is 
selected). A global solar radiation frequency histogram is obtained 
by dividing the number of observations in each bin by the total 
number of observations in the data set. As shown from Fig. 6, the 
highest peak frequencies corresponds to the range of global solar 
radiation between 0 and 50 W/m? with a value equal to 54.15%. 
Concerning the other ranges of solar radiation, the highest 


frequency is related to the bin 250-300 W/m? with a value of 
3.39% while the lowest one corresponds to the range of solar 
radiation within 950-1000 W/m? with a value equal to 0.11%. We 
can observe that all the other range of solar radiation have near 
values of frequency distribution. 

Fig. 7(a)-(d) show the seasonal frequency distributions of 
global solar radiation based on the data of the three last years 
(2008, 2009 and 2010). Concerning the winter season, according 
to Fig. 7(a) and ignoring the solar radiation within 0 and 50 W/m? 
which corresponds almost to the night periods, it is observed that 
the prevalent range of global solar radiation is related to the range 
of 200-250 W/m? with a probability of 5.23%. It is worth noting 
also that during this season the radiations greater than 550 W/m? 
have a very low occurrence probability due to the cloudy periods 
of this season. Fig. 7(b) shows the frequency distribution of the 
global solar radiation in the summer season. It is observed that 
the highest frequency is related to the range of solar radiation 
750-800 W/m? with a value of 4.48% while the lowest frequency 
corresponds to the radiation values between 950 and 1000 W/m7. 
It is noticeable that all ranges of global solar radiation have 
significant values of occurrence probability due to the significant 
sunny months of this season. As shown in Fig. 7(c), the frequency 
of the global solar radiation during the autumnal months follow a 
similar distribution pattern but the peak of frequency corre- 
sponds to the solar radiation which falls in the range of 50- 
100 W/m? with a value equal to 4.21% because this radiation can 
be registered both during the sunrise and the sunset while the 
important radiations such those within 950 and 1000 W/m? can 
be recorded only during the hours surrounding the noon that is 
why their occurrence probability are very low compared rela- 
tively to the lowest radiations. During this season, the registered 
global solar radiation did not exceed the value of 800 W/m? and it 
is attributed mainly to the cloudy sky conditions of the autumnal 
months and the elevation of the sun. Fig. 7(d) shows that during 
the spring season the global solar radiation spread reaching the 
1000 W/m? in comparison relatively to the autumn and winter 
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Fig. 7. Frequency distribution of global solar radiation during winter season (a), summer season (b), autumn season (c) and spring season (d). 


seasons. On can observe also that with ignoring the solar of 
radiation that falls in the range of 0-50 W/m’, the prevalent solar 
radiation corresponds to whose falls in the range of 550-600 W/ 
m? with a probability of 4.08% while the infrequent solar radia- 
tions are related to whose falls within 900 and 1000 W/m? with a 
frequency of 0.84%. The obtained outcomes give good information 
about the solar radiation availability in the studied site and 
confirm again that the gulf of Tunis presents an important solar 
potential. 


5.3. Comparison of results and test of model 


Fig. 8(a)-(d) show the comparison between the measured 
global solar flux and the predicted one using a conventional 
model during four specific days. We note that these specific days 
correspond to the summer and winter solstice days and the vernal 
and autumnal equinox days. It is found that the conventional 
model consistently over-predicts the measurements at all times. 
This is expected since the model does not account for local 
weather conditions such as the presence of cloud and dust. 
Moreover, one can observe that there are visible discrepancies 
between the measurements for the three years. The performance 
of the conventional model was evaluated in terms of the relative 
percentage error (MAPE).Table 3 summarizes the MAPE of each 
year (2008, 2009 and 2010) as well as the global error for the 
whole samples during the four above-mentioned days. It is found 
that the highest global MAPE is occurred in the autumnal 
equinox, winter solstice and vernal equinox days with values of 
14.26%, 14% and 10.34% respectively. This is generally attributed 
to the local cloud formations in these wintry, autumnally 
and vernally days which obviously vary from year-to-year. 
In contrast, the results presented for the summer solstice day in 
Fig. 8(b) show that both the measurements and the conventional 


model predictions differ only slightly. In fact, there is a practically 
close agreement between the results of the conventional model 
and the mean values of the measurements with a slight over- 
prediction of the global solar radiation received by a horizontal 
plane by the conventional model with a global MAPE equal to 
4.10%. Results for the other days of the years, which are not 
presented here follow a similar trend and agree with the above 
conclusion, i.e. some discrepancies occur among the measure- 
ments for different years and with the predictions for the cloudy 
periods, namely, October-May. These differences decrease 
obviously for the months of June-September (summer months) 
for which the sky is relatively clearer. 

The promising similarities between the measured and pre- 
dicted solar-flux variations during the clear-sky dates suggest 
that global solar radiations have to estimated basing on this 
conventional model but it is will be more interesting to focused 
on refining a specific theoretical model concerning the gulf 
of Tunis. 


6. Conclusion 


The determination of global solar radiation as well as its real 
distribution at any site is vital important for many scientific, 
engineering and environmental applications. In this context, this 
paper is subscribed by the assessment of the solar energy 
potential in the gulf of Tunis, Tunisia. In fact, the ambient 
temperature, the global solar flux, the sun duration, the yearly 
and the seasonal frequency distribution of the global solar 
radiation were presented. The results proved that the site of 
Borj-cedria in the gulf of Tunis has been endowed by a promising 
solar energy potential. A conventional model of the global solar 
flux estimation was validated by the measurements data in 
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Fig. 8. Comparison between the predicted and the measured Global solar radiation during the winter solstice (a), summer solstice (b), vernal equinox (c) and autumnal 


equinox (d) days. 


Table 3 
Relative error percentage between the daily predicted and the measured global 
solar radiation. 


REP 2008 REP 2009 REP 10 Global REP 
(%) (%) (%) (%) 
Winter solstice day 15.89 18.14 8.56 14.2 
Summer solstice day 3.49 2.66 6.16 4.1 
Autumnal equinox 14.49 12.79 15.49 14.26 
day 
Vernal equinox day 11.36 6.18 13.75 10.34 


critical days of the last three years. A daily comparison basing on 
the mean absolute percentage error showed that this model 
predicted with good accuracy the real data under the clear sky 
conditions with an error of 4.10% while its performance decreases 
in the cloudy conditions with an error reaching the value of 
14.26%. Hence, it is recommended to develop an appropriate site 
model in order to predict satisfactory the global solar flux in the 
gulf of Tunis. Finally, this aim is very important because if there 
are missing data for some months due to faulty operation of the 
measuring device, the developed model would be very useful 


used to predict the missing monthly global solar radiation at any 
date in the gulf of Tunis. 
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